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ABSTRACT 



To reconstruct dark energy models the redshift Zgq, marking the end of radiation 
era and the beginning of matter-dominated era, can play a role as important as Zt, the 
redshift at which deceleration parameter experiences a signature flip. To implement 
the idea we propose a variable equation of state for matter that can bring a smooth 
transition from radiation to matter-dominated era in a single model. A popular A oc p 
dark energy model is chosen for demonstration but found to be unacceptable. An 
alternative A oc pa^ model is proposed and found to be more close to observation. 

Subject headings: Cosmology; Variable A 
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Introduction 



The current trend in modern cosmology is mainly focussed on issues of possible explan ations 



of the observed 



Perlmutter et al. 



ate t ime acceleration in the expansion of our Universe (iRiess et al. 



1998 



1999I ). The presence of a dark energy of unknown origin which dominates over 
matter in our Universe is now almost accepted. Unlike ordinary matter the dark energy has 
a repulsive effect and its dominance accounts for the observed acceleration. Several candidates 
of dark ene rgy are found in lite r ature. Among them the most elementary is t he cosmologica 



constant AdPadmana 



1992 



2004 



Wetterich 



1995 



lan 



2003 



Arbab 



Dutta Choudhury fc Sil 



Peebles and Ratra 



997 



Padmanabhan 



2001 



20031). Variable A mode l sdCarvalho et a. 



Vishwakarma 



2002 



Shapiro fc Sola 



20061 ) were introduced later in an effort to overcome the fine 



tuning probl em. Among ot 



ler c andidates m ost popular are d ifferent kinds of sca lar fields like 



quintessence (Icaldwell et al.l 



19981 ). K-essence( 



to name only a few. The Chaplygin gas llKamenshchik et al. 



Chiba et al. 



20001) and phantom fields (Icaldwell 



2002) 



20011 ) has an equation of state that 



can produce negative pressure in later phase of evolution to account for the accelerated expansion 
of the Universe . Some geometrical origin of the late-time acceleration are also found in modified 



gravity models llNojiri and Odintsov 



20071 ). All these models have their own merits and demerits. 



The exact nature of dark energy is therefore yet to be a scertained. 



different dark energy candidates one may see the work of 



Sahnil fl2004h or 



^or a detailed review on 



Copeland et al. 



fl2006h . 



Most of these models are usually studied only in the matter dominated phase assuming the 
presence of some kind of dark energy to explain the transition from deceleration to acceleration. 
This is due to the fact that radiation era was very short lived and dominated the Universe in its 
early phase of evolution whereas acceleration is a more recent phenomena which has occurred 
very late in the matter dominated epoch. The viability of such models are usually checked by 
comparing the variation of deceleration parameter q with the red-shift z obtained theoretically 
with that from observations. A good estimation oi z = Zt, the redshift marking the transition 
from deceleration to acceleration, is available from observation. But observational data are 
limited over a small range of z values around Zt. Hence many models that fit observational data 
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are found to behave quite differently outside this range. To choose among these models we need 
at least another event that happens to be outside this range of 2;, preferably in the distant past. 
We point out that a very good estimate of z = z^q, the redshift when matter density was equal 
to radiation density, is at our disposal. Thus instead of a single point we now have two well 
separated points through which the q vs. z curve for any viable cosmological model should pass. 
Such a model would be more reliable as its viability could be ensured over a large range of z 
values. But then to estimate z^q we must have both radiation and matter in a single model. 

In the next section we propose a variable equation of state for matter that mimics radiation 
in the past and pressureless dust in the later course of evolution. In section 3 we examine the 
A oc p model with the proposed equation of state for matter. In section 4 an alternative variable 
A model is proposed. 



2. Radiation and matter in a single model 



To begin with let us keep aside dark energy for a while and follow the standard practice of 
idealizing the Universe as homogeneous and isotropic supported by the isotropy of microwave 



background radiation. The spatial ge ometry of the Ro 



view of strong observational evidences (I Zhao et al 



Dertson- Walker spacetime is chosen flat in 



20071 ). In these kinds of standard models of the 



Universe with no dark-energy, rad iation domina tes for roughly the first 2000 years of evolution 



and matter dominates there- after (jPeacock 



2OO3I ). In Friedmann models the homogeneous and 



isotropic Universe is filled with an ideal fluid. The equation of state for this fluid is either that 
of radiation {pr = |pr) or that of pressureless dust {pm = 0). Thus one needs two different 
cosmological models with different equations of state for matter to explain the expansion history 
in two different epochs of the same Universe. No single Friedmann model can describe the full 
history of the Universe. The deceleration parameter q also takes different positive values, however 
constant, in each epoch. Such cosmological models are usually studied by joining a radiation 
Universe (j9 = p/3) to a matter dominated Universe (p = 0), smoothly at some particular time. 
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A ttempts w e re made to describe the entire evoluti o n in a s i ngle m odel by IChernin 



19681), 



McintoshI (11967 . 



19681) 



Landsberg fc ParkI (119751 ) 



Jacobs! (119671 ) and 



Cohenl (119671 ) in late 



(119651 . 



sixties of the last century. In those models the homogeneous and isotropic Universe is assumed 
to be filled with a perfect fiuid which has two components viz. radiation and matter. Assuming 
no interaction between radiation and matter, the radiation energy density pr decreases with time 
in the usual manner ~ while the matter density falls as pm ~ cl~^- Thus the total energy 
density p and the total pressure p of the fiuid may be written as 



P = Pr + Pm = Proiao/a)^ + Pmoiao/af 

P = Pr+Pm = ^Pro(ao/a)^ 



(1) 

(2) 



where the subscript '0' indicates the value at present epoch. We thus find that p and p maintain 
an equation of state 



p 1 

/ = - = o [1 + iPmo/Pro){a/ao)] 
P 3 



-1 



1 + ^ 



where Ggq is defined as the scale when pr and pm are equal i.e. 

Proi'^O / (^eqY = Pmoi'^0 / (^eq)^ ■ 



(3) 



(4) 



One may see that for small values of 'a' the equation of state approximates to that of pure 
radiation while for a >> a^q the pressure becomes negligible and mimics the matter dominated 
epoch. 

We propose the cosmological fiuid to obey the above equation of state. With such a model 
in presence of dark energy it becomes possible and easier to tune parameters of the model in 
order to obtain desired or observed values for both z^q and zt simultaneously. As an exercise, in 
the next section, we shall consider a variable A 'dark energy' model to demonstrate this point. 



3. Testing a A oc p model 

We need an example of a dark energy model to show how our idea can be implemented. 
In absence of any knowledge about the exact nature of dark energy majority of approaches are 
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phenomenological where we limit the equation of state for dark energy from observational data. 
Althou gh observation does not rule gut the possibility of a varying equation of st ate for dark 



energy (IKujat et al 



2002 



Bartelmann 



2005 



Knop et al 



2003 



Tegmark et al 



20041 ) the evidence s 



are s trongly in favour of an equation of state that is very close to PdeI Pde = — KiDavis et al. 



20071 ). Thus observational data supports the choice of Einstein's cosmological constant A as 
a dark energy. However, a major challenge in constant A cosmology is perhaps to explain its 
vanishingly small value at present. From field theoretical point of view if A repr esents the vacuu m 



19891) 



energy density then its value should be ~ 10^^^ times its observed value today (iWeinberg 
Attempts were made to explain such anomalies in the order of magnitude of A by introducing 
the idea of a time varying A. 



There exist 


s a large num 


3er of variable A models in the literature 


fCarvalho et al. 


1992 


Wetterich 


1995; 


Arbab 


1997; 


Padmanabhan 


2001 




Vishwakarma 


2002 


. 


Shapiro & Sola 


2004; 


Dutta Choudhurv & Sil 


2006 


Rav. Mukhopadhvav and Dutta Choudhurv 


2007 


). A popular 



choice for variable A in the literature is a function that either decays like the energy density p 
or like the squar e of the Hubble parameter or as d/a . However, all these models were found 



to be equivalent (IRay. Mukhopadhyay and Meng 



20071 ). A difficulty in such models is to find 



a 'single' positive constant of proportionality that can produce deceleration in the radiation 
dominated era as well as a transition to acceleration late in the matter dominated era. This is 
indeed obvious (though undesirable) as in such models the variations of the matter density p are 
distinctly different in radiation and matter dominated eras. Thus with a variable equation of 
state for matter that smoothly connects the radiation and matter era in a single model, it might 
be interesting to estimate the value of the constant of proportionality, that produces acceptable 
values for observational parameters. In this section we choose a variable A oc p as the dark energy 
candidate and investigate how far it can produce observational results including a desired value 
of Zeq when we introduce a variable equation of state for matter. Keeping in mind the above 
said difficulty we will relax the reasoning behind the choice of the equation of state for matter 
by introducing a parameter m in the expression for / in equation ([3]) preserving its asymptotic 
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behavior 



/ 



1 + 



"eq 



m-\ —1 



(5) 



where m is a positive constant that may be determined from observational constraints. 



With a choice of the energy-momentum tensor of the Universe 



8^ 



g^"' = {p + p)u^u'' - pg 



c^A 
8^' 



flU 



we write the Einstein's field equations 



2 ^ 



(6) 



(7) 



for a homogeneous isotropic Robertson- Walker type spacetime with flat {k = 0) spatial geometry 
as 



a 
a 



— = 2up+ -c'A 
d o 

-u{p + 3p) + ^c^A 



(9) 



where u = ^p. 



From equation ([S]) it may appear that with A oc p, one can trivially absorb A in to the 
constant u and hence there would be no dark energy in the cosmological evolution. This is indeed 
true for a perfect fluid with an equation of state p = up. The introduction of A oc p effectively 
modifies the proportionality constant n and the model becomes equivalent to one with a new 
fluid with no dark energy. But with a choice of variable equation of state for matter one does 
not have the same liberty in equation (Q. Following Bianchi identity, the divergence of the field 
equation ([7]) gives, 



[R 



-Rg 







-SttG 



fJ,U 



c^A(t) 
SttG 



9'% 



(10) 



which indicates that empty spacetime cannot be obtained as solutions of the field equations if 
A varies with time. Thus a variable A is consistent w ith Mach's principle that for a meaningful 



spacetime geometry presence of matter is a necessity ljVishwakarma 
energy-momentum tensor requires 

up + ^Ac^ + 3m(p + p)- = 0. 
D a 



2002). The conservation of 
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The above equation clearly shows the presence of interaction between matter-radiation and dark 
energy. 

Among the set of three equations ([8]), ([9]) and (fTT!) only two are independent, whereas 
number of unknown functions are four. In order to get an exact solution to the problem we 
make the following choices. We make an ansatz that the effective equation of state for the 
matter-radiation part as given in equation ([5]) is 



/ 



1 + 



eq 



and choose the variable cosmological parameter as 

u 



A 



;Xp. 



'12) 



Here A is a positive constant. For the above choices equations 

d2 A~ 



a 
a 



and (IHl) reduce to 



— = up{2 + - 



The deceleration parameter q is therefore given by 

a/ a 1 + 3/ - I 



jo? 



2 + 1 



(13) 
(14) 

(15) 



Choice of / as in equation ^ implies that for small values of a the deceleration parameter 
behaves as If we believe that in early radiation dominated era the Universe went through 
a decelerated phase with g ~ 1 then A must be vanishingly small in the radiation era. In the 
later phase for a >> Oeg, the function / becomes very small and q behaves as tt-t- Hence to 
observe late-time acceleration in the expansion we must have A > 3, preferably A > 12 so that 
q < —0.5 for large a to make it consistent with observation. Restriction on A also comes from the 
observation of density parameters Qm and Q\ defined as 

2up 6 



nr. 



c'A 



uXp 



6 + A 
A 



3dVa2 3dVa2 6 + A' 



(16) 
(17) 
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Obse rved values of Qm and at the present epoch are Qmo ~ 0.3 and Q^o ~ 0.7 flRiess et al 



20011 ). Such observation puts a hmit on the present value of A as Aq ~ 14. 



With the expression for / as in equation ([5]) the expression for deceleration parameter q in 
view of equation (ITSl) becomes 

3[(f-r + l]-i + (3-A) 



6 + A 



18) 



It is evident from the above equation that to ensure a decelerated phase in the beginning of 
evolution and a late time transition from decelerated to accelerated phase in the expansion of 
the Universe we require 3 < A < 6. If the value of A lies within this range one can see that the 
range for q will be — | < g < |. We can see that this makes it difficult to produce observational 
requirements. The radiation-matter density p can be solved analytically from conservation 
equation (fTTl) to get 

p = A[a-^(a- + ar,)^]^ (19) 

where A is the constant of integration. Notice that for A = (i.e no A parameter) and for any 
value of m > the energy density p ~ for small values of a (such that a << a^q) whereas for 
large values of a (such that a >> a^q), p ~ a~^. Substituting (|T9l) back in ([8]) we can solve for a 
for different possible values of m. For a choice, m = the equation can be easily solved to 
obtain 

a=[{at + af )l -a":^]^ (20) 
where, a is a constant involving A and A 



9 



a 



(A + 6) 



Au{2 + ^ 



(21) 



Note that the particular choice m 
from equation (l20l) we get 



12 



(A+6) 



gives m = 2 for A = (i.e no A parameter) and hence 



~eq 



at 

1 + — 



— a 



Cteq , 



eq 



(22) 



For such a model with no A parameter the scale factor a ~ for small values of t (radiation 

2 

dominated era), and for large values of t (matter dominated era), a ^ t^. Thus instead of two 
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different (radiation dominated and matter dominated FRW) models we can have a single model 
with variable equation of state that h as the desirable asymptotic behavior. This finding may be 
identified with those of I Jacobs! (Il967l ). 



It may be noted at this point that in general there exists no relation between m and A. If 
we do not go for any particular choice of m it might be interesting to see how the observational 
constraints can relate them. If a = at marks the transition from deceleration to acceleration in 
the expansion of the Universe we can see from equation f|T8|) that 



at 



6- A 



(23) 



A-3 

Since the acceleration of the Universe has occured late in the matter dominated era, the above 



ratio must be gre ater t 



lan 1 and hence A < 4.5. This resets the limits on A as 3 < A < 4.5. 



Riess et al 



(j200l|) and 



Turner &: Riesd (120021 ) have shown the redshift- distance relation of SN 
la compared to a family of flat, cosmologies where the transition between accelerating and 
decelerating epochs occurs at a redshift 



Zt 



2fi; 



1. 



For ~ 0.7 and 0.^ ~ 0.3, the above relation gives 

ao 



l+Zt = 



1.7. 



at 



Al so the redshift Zf.q marking the beginning of matter dominance over radiation is usually given 



by ([Peacock 



20031 ) 



1 + 2;, 



eq 



On 

— ^ 1000. 

Oieq 



Assuming ^ ~ > '^^ fi^d that m and A obey the relation 

■}^y\>i^, (24) 

1.7/ A-3 ^ ^ 

We use equation (l2ll) to eliminate A from equation (fTSll . The present day value of the deceleration 
parameter go is now a function of m as given below 



?0 



103m)(4 + 3(^)3m)- 



(25) 
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When we plot the above function in fig.([T]) we find there exists a minima at m ^ 0.242 (and hence 
A ~ 3.528). The minimum value of go is however very small (~ —5.623 x 10^^). There exists 
no other choice of the pair m and A for which go could be made smaller than this value. Thus 
we find that a constant A that fits both the radiation and matter dominated models is capable 
of producing a signature fiip in the deceleration parameter. However within the framework 
of our model the value of go becomes insignificantly small due to observational constraints 
on the values of and a^. The variation of g with a is plotted in fig. ([2]) in ao = 1 unit for 
a pair of m and A that produces maximum possible acceleration at the present epoch in this model. 



- 0.001 

- 0.002 
3 - 0.003 

- 0.004 

- 0.005 



Fig. 1. — go vs TO plot. The minima in the graph shows that go can never be made smaller than —5.623 x 10'^. 

0.1 
0.08 
0.06 
0.04 
0.02 


- 0.02 

0.2 0.4 0.6 0.8 1 

a 

Fig. 2. — q vs a plot. The plot shows the variation of q with a for to = 0.24207 and A = 3.528. 
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4. An alternative variable A model that passes the test 

In view of the failure of a variable A model of the kind A oc p to meet observational 
constraints as well as desired features as shown in the previous section we propose that a simple 
choice for variable A satisfying the observational constraints might be 

A = ^A(a)p (26) 

where 

A(a) = 3(^) ■ ^2^) 

The motivation behind suc h a choice is the success of 'Lambda Cold Dark Matter' (ACDM) 



mo del (IPadmanabhan 



20031 ) in explaining the observed acceleration. Anticipating that in the 
matter dominated era the energy density p would vary as a~^, we propose fl26l) so that A may 
behave as a cosmological constant in the later phase of evolution. This way the field theoretical 
demand of very large value of A in the early Universe will also be met if p ~ for very small 
values of a. 

Using equation ( 1271) and equation ([3]) the deceleration parameter q, given by equation (|T5l) . 
in this model takes up the following functional form 

1 + —^ ' + 3 



<1=- 7^T^-1- (28) 

The asymptotic behavior of q shows that it varies smoothly from +1 in distant past to —1 in 
distant future with its signature flip occuring at a ~ at- When q is plotted against a (oq = 1 unit) 
in fig. ([2]) we find that for a considerable period of time in the early stage of matter dominated era 
q remains almost constant at a value +0.5. This is indeed a very desirable feature of this model 
as it provides a favourable condition for structure formation. A rough estimation with Zgq ^ 1000 
and Zt ^ 0.7 shows that the present value of deceleration parameter go ~ —0.56. 



The conservation equation (fTTj) with the above choice of A reduces to 



P 
P 



3A 
6 + A 



+ 2(g+l) 



^ = 
a 



(29) 
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Fig. 3. — q vs a plot. 



which imphes 




In p = In A + In a'^ + In (a^g + a)^" + In (2af + a^)'^'^ + In - + a 



2N27 



- In (2^/^ + a)^^ + r] arctan 



2a 



(30) 



where A is the arbitrary constant of integration and a, (3, 7 and rj are constants involving at 
and aeq. For Oeg = 0.001 and at = 0.59 in oq = 1 unit, the values of the constants are a = 0.5, 



P = 4.06 X 10 



-10 



7 



1.12 X 10 ^ and 77 = 7.76 x 10 ^. Thus only the second and the third 



term on the right hand side of equation fl30l) contribute significantly to the variation of p with a. 
Variation of p and A with a is shown below [fig.(jl])] in arbitrary units. The density parameters 
are also plotted in fig.®. In both the figures a is scaled in = 1 unit. 



The asymptotic behavior of p may be easily investigated from equation (!29|) . For very small 
values of a we find that A ^ while g — > 1. This makes the conservation equation (1291) to reduce 
to 

^ + 4^ = 0. (31) 
p a 

Hence the matter-radiation density p ~ in the early phase of evolution whereas for a — > cxd 
the functions A ^ 00 while g — — 1 making equation (12^ to reduce to 

pa , ^ 

- + 3- = 0. (32) 
p a 
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Thus in the late stage of evolution the matter-radiation density p a ^. This makes A to attain 
a constant value in the later course of evolution. 




0.2 0.4 0.6 0.8 1 



a 

Fig. 4. — p , A vs a plot. In the early stage of evolution A ~ as p ~ a~^. A gradually attains a constant 
value with evolution since p for large a. 




a 



Fig. 5. — Qm , vs a plot. 

5. Conclusion 

We have shown in this present article that the viability of a dark energy model can be 
cheked in a better way by demanding that Zeg should be ~ 10^. Introducing a variable equation 
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of state for matter we have investigated the nature of expansion of the Universe in two variable 
dark energy models. The variable equation of state for matter ensures that the early epoch of 
the Universe is dominated by radiation. With the expansion of the Universe the pressure falls as 
a~"^ (m > 0) to produce an almost pressure-free late-Universe. The equation of state for dark 
energy is essentially chosen to be aoE — Pde/ Pde — —1 while the poE is in general a function 
of the scale factor a. This in other way amounts to the same as the presence of a dynamical A 
in the field equations. It makes the conservation equation notably different. One has to allow 
interaction among the matter part and dark energy. The functional form of A is then chosen 
to produce late-time acceleration in the expansion of the Universe. A A oc p model fails to 
meet observational constraints. A choice of A cx pa? is found to be a better alternative. Thus 
we conclude that meeting the additional Zeq ~ 10^ requirement may be an essential feature of 
reconstructing dark energy models. 
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